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IS.  Abstract 

The  attenuation  of  an  acoustic  disturbance  propagating  In  a  medium  which 
is  not  perfectly  elastic,  homogeneous  and  iaotroplc  can  be  considered  to  be 
the  result  of  energy  dissipation  due  to  classical  processes  and  scattering. 

The  contribution  of  each  of  these  mechanisms  Is  dependent  upon  the  frequency 
of  the  sound  wave.  In  this  report  several  mathematical  models  are  discussed 
for  a  transverse  wave  moving  In  a  visco-elastic  material  and  a  scattering 
theory  Is  developed.  The  theory  is  applied  to  sea  ice  and  attenuation  as  a 
function  of  frequency  Is  predicted  for  this  material. 


NATIONAI  If'MNICA! 
INFORMATION  TPV'i  F 

u  s  rw*prt*tr,-f”*  •  <  * 

Spring,** i  i  .  .  IM 


IT.  Key  Ward* 

Acoustic  attenuation  In  sea  ice 
Sea  Ice,  acoustic  attenuation  in 
Transverse  acoustic  waves 


IS.  DleWibaMaai  ttatom *w> 

Document  is  available  to  the  public 
through  the  National  Technical 
Information  Service,  Springfield, 
Virginia  22151 


It.  Security  Clattif.  (•<  Silt  mpartT 

Unclassified 


S.  Saaactty  Claetil.  (W  Ala 

Unclassified 


21.  Ha.  ,1  Pagat 


n.  Net 


SgptadwcHow  *♦  oaamlatad  papa 

S 

/ 


Pona  DOT  P  1700.7  (S-TO 


autfcaHgad 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


CONTENTS 

INTRODUCTION .  1 

LONG  WAVELENGTH  ISOTROPIC  SOLUTIONS .  2 

Derivation  of  the  General  Equations .  2 

Mathematical  Models  of  Visco-Elastic  Materials .  1 

Maxwell  Model .  4 

Voigt  Model .  3 

Expansion  to  More  Complex  Models . .  .  10 

Extension  of  the  Theory  to  Include  Anisotropy .  15 

SHORT  WAVELENGTH  SOLUTIONS .  16 

General  Discussion .  16 

Scattering  of  a  Plane  Wave  from  a  Spherical  Obstacle .  17 

APPLICATION  OF  THE  THEORIES  TO  SEA  ICE .  23 

Scattering  Attenuation  Approximation .  23 

Total  Attenuation  Approximation .  28 

Depth  Dependence  of  Shear  Wave  Propagation .  29 

CONCLUSIONS .  32 

APPENDIX  A,  Description  of  In-Situ  Shear  Wave  Experiment .  35 

APPENDIX  B,  Computer  Programs . 37 

REFERENCES .  42 


APL-UW  7310  iii 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


INTRODUCTION 


Activities  in  much  of  the  arctic  region  require  a  knowledge  of  ice  thick¬ 
ness  to  assure  safety  of  equipment  and  personnel.  This  knowledge  is  especially 
essential  during  exercises  that  require  the  operation  of  heavy  equipment,  for 
example,  the  employment  of  aircraft  on  ice  runways  for  logistical  support. 

During  the  last  75  years,  many  studies  of  ice  have  been  performed  to  under¬ 
stand  its  physical ,  mechanical  and  chemical  properties.  A  thorough  discussion 
of  this  work  has  oeen  published  by  Weeks  and  Assur1,  so  an  extensive  literature 
review  will  not  b3  presented  here.  Even  though  some  of  these  studies  show  ice 
to  be  a  relatively  good  conductor  of  acoustic  disturbances,  an  efficient  method 
of  acoustically  determining  ice  thickness  from  the  top  surface  still  does  not 
exist.  Instead,  ice  thickness  measurements  are  presently  based  on  other  physi¬ 
cal  characteristics  of  ice  and  water,  such  as  dielectric  permittivity,  dielec¬ 
tric  loss  tangent  in  ice  and  in  sea  water,  or  elastic  oscillations  of  the  total 
ice  cover. 

Simple  calculations  based  on  the  idealized  assumption  of  isotropic,  homo¬ 
geneous  ice  indicate  that,  because  of  acoustic  impedence  matching  at  the  ice- 
water  interface,  au  incident  acoustic  compressional  vive  will  be  almost  entirely 
transmitted  into  the  water  medium.  Conversely,  because  the  water  will  not 
support  shear  waves,  the  shear  wave  stress  is  assumed  to  vanish  at  the  boundary. 
This  assumption  indicates  that  an  acoustic  shear  wave  impinging  on  the  boundary 
will  be  reflected  totally  at  normal  incidence  and  less  at  greater  angles.  Thus, 
neglecting  real  and  apparent  attenuation,  it  would  appear  from  this  simple  model 
that  the  best  method  of  determining  the  ice  thickness  acoustically  from  the  top 
surface  would  be  to  generate  a  shear  wave  in  the  ice  and,  knowing  the  wave  speed, 
calculate  the  thickness  from  the  reflection  of  the  wave  from  the  ice-water 
boundary . 

It  is  known,  however,  that  sea-ice  is  not  a  simple  medium,  but  rather  a 
highly  complex  inhomogeneous,  anisotropic,  polycrystalline  material  for  which 
the  acoustic  propagation  characteristics  may  be  dependent  on  the  growth  and 
life  history.  Therefore,  the  ideal  model  discussed  above  is  a  gross  oversim¬ 
plification  of  the  physics  of  the  problem.  One  principal  area  of  concern  that 
must  be  resolved  before  ice  thickness  can  be  reliably  measured  with  instrumen¬ 
tation  based  on  shear  wave  propagation  is  the  attenuation  of  the  wave  in  the 
medium. 

The  attenuation  of  a  sound  wave  propagating  in  a  material  that  is  not  per¬ 
fectly  elastic,  homogeneous  and  isotropic  can  be  considered  to  be  the  result  of 
two  mechanisms  which,  for  the  purposes  of  this  report,  will  be  referred  to  as 
classical  dissipation  and  scattering.  Classical  dissipation  originates  from 
many  sources,  a  few  of  which  are  internal  friction,  anelastic  behavior  of 
the  material,  thermal  dissipation,  viscous  slippage  at  crystal  boundaries  and 
the  presence  of  dislocations  and  solute  atoms  in  the  medium.  Energy  losses 
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due  to  scattering  originate  from  the  interaction  of  the  acoustic  wave  with 
scattering  centers  in  the  medium. 

The  relative  contribution  of  each  of  these  mechanisms  to  the  total  attenu¬ 
ation  depends  c&  the  frequency  of  the  sound  wave.  At  low  frequencies,  the  wave¬ 
length  of  the  sound  is  very  large  compared  to  cne  scattering  centers;  therefore, 
the  scattering  cross  section,  i.e.,  the  relative  amount  of  energy  scattered  out 
of  the  incident  wave,  is  extremely  small  and  the  attenuation  is  due  almost 
entirely  to  classical  dissipation.  As  the  frequency  increases,  the  attenuation 
due  to  scattering  becomes  more  important  until,  when  the  wavelength  becomes 
approximately  on  the  order  of  the  size  of  the  scattering  center,  scattering 
eventually  predominates. 

We  have  performed  a  calculation  of  the  attenuation  from  both  classical 
dissipation  and  scattering  which  ultimately  will  allow  the  prediction  of  the 
total  attenuation  of  a  shear  wave  in  sea  ice.  This  report  gives  the  results 
of  that  calculation  and  recommends  that  experiments  be  performed  to  evaluate 
the  theories. 


LONG  WAVELENGTH  ISOTROPIC  SOLUTIONS 

DERIVATION  OF  THE  GENERAL  EQUATIONS 

For  wavelengths  that  are  long  compared  to  the  size  of  the  scattering 
centers,  i.e.,  at  low  frequencies,  the  calculation  of  acoustic  velocity  and 
attenuation  can  be  derived  ignoring  the  contributions  due  to  scattering.  Al¬ 
though  all  of  the  internal  mechanisms  contributing  to  classical  dissipation 
are  unknown,  the  general  theory  can  be  developed  by  grouping  all  attenuation 
sources  as  a  simple  constant.  This  is  accomplished  mathematically  by  the  re¬ 
placement  of  the  shear  modulus  (M)  and  the  Lame'  constant  (A)  in  the  elastic 
stress-strain  relationships  by  the  first  order  differential  operators 

Ms  w'k 

3  '  il] 

A'  X  + 

In  Eq.  (1),  the  unprimed  terms  denote  the  elastic  and  the  primed  terms  denote 
the  viscous  (or  attenuative)  contributions. 
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Consider  an  acoustic  shear  wave*  moving  in  the  positive  (downward)  x- 
direction.  The  displacement  potential  is  given  by 


♦  ,  AeU“t-ax) 


(2) 


where 


u)  is  the  circular  frequency 

t  is  the  time 

a  is  the  complex  wave  number. 

This  displacement  potential  obeys  the  wave  equation 

P  =  M  V2  \p  ,  (5) 

3t 

where  p  is  the  density  of  the  medium  and  V2  is  the  Laplacian  operator. 
Substitution  of  Fqs.  (1)  and  (2)  into  Eq.  (3)  yields 

p  =  yVVy'V2  (4) 

dt2  3t 

or,  carrying  out  the  operations, 

pu)2  =  a2  (u  ♦  iu>p')  .  (5) 

For  a  dissipative  medium,  the  complex  wave  number  can  now  be  written  as 

a  =  k  -  it  ,  (6) 


where 


k  is  the  running  vector,  u/c 

c  is  the  wave  velocity 

t  is  the  generalized  attenuation. 


Thus,  writing  the  potential  in  the  form  of  Eq.  (2), 


Ae-Zxeiu(t  '  x/c) 


(7) 


*The  derivation  for  the  compressional  wave  propagation  is  similar  to  that 
presented  here,  the  difference  being  the  introduction  of  the  (A  ♦  2M) 
operator  in  the  wave  equation,  Eq.  (3),  rather  than  the  shear  modulus  (M) 
operator.  The  displacement  potential  must,  of  course,  be  associated  with 
the  compressional  rather  than  the  transverse  wave. 
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which  denotes  a  wave  propagating  in  the  positive  x-direction  with  a  velocity 
c  and  an  amplitude  which  decays  at  a  rate  of  t  per  unit  distance  as  the  wave 
travels  into  the  medium. 


Making  this  substitution  into  Eq.  (5)  and  separating  the  real  and  imaginary 
components  yield  the  simultaneous  equations 

pu)2  =  u(k2  -  i2)  +  2p "coki 

0  =  wi/(k2  -  i2)  -  2pkT  .  (8) 


Setting 


these  equations  have  the  solutions 


k2 


J_  pu)2 
2  ~ 


[>/i  ♦  rV  ♦  i 

1  l  ♦  rV 


(9) 


where  now  both  p  and  R  may  be  frequency  dependent.  The  procedure  from  this 
point  is  to  introduce  a  mathematical  model  ot  the  medium  composed  of  elastic 
and  viscous  constants  which  are  frequency  independent,  and  then  fit  the 
results  to  experimental  data. 


MATHEMATICAL  MODELS  OF  VISCO-ELASTIC  MATERIALS 

Historically,  there  are  many  different  rheological  models  that  have  been 
used  to  represent  viscous  materials.  Although  these  models  have  been  devel¬ 
oped  to  interpret  the  results  of  static  measurements,  one  or  more  of  them  can 
be  applied  to  the  results  of  dynamic  tests,  particularly  if  limited  frequency 
ranges  are  being  considered. 

Maxwell  Model 


Consider  the  mechanical  representation  of  a  visco-elastic  medium  shown  in 
Figure  1.  Here,  the  material  is  represented  by  a  spring  denoting  the  elastic 
element  (Em)  in  series  with  a  dash-pot  denoting  the  viscous,  or  dissipative, 
element  (hm)  •  In  this  model  both  Era  and  Pm  are  frequency  independent. 
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Figure  1.  Maxuell's  Mechanical 
Representation  of 
Visco-Elastic  Solids. 


Maxwell2  defined  the  stress-strain  relationship  for  this  visco-elastic 
solid  as: 


.  .  a. . 

o.  .  =  2M  y-  S. .  --U  (i  +  j) 

dt  it  dt  it  5  v  r 


(10) 


where  Ojj  and  Sjj  are  components  of  the  stress  tensor  and  the  strain  tensor, 
respectively,  and  6  is  the  relaxation  time  of  the  material.  If  the  wave  equa¬ 
tion  is  to  be  applicable,  Hooke's  law  states  that  the  relationship  between 
stress  and  strain  must  be  of  the  form 

o  =  MS  ,  (11) 

which,  for  this  case,  may  be  accomplished  by  writing  the  shear  modulus  (M)  as 
an  operator, 


M  = 


y 


(12) 


or,  in  terms  of  the  elastic  and  viscous  constants  of  Figure  1, 


M  = 


(13) 


Thus,  substituting  Eq.  (13)  into  E  (11)  and  simplifying  yield  the  relation¬ 
ship 


d 

dt 


m 


E  S 
m  dt  m 


(14) 
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Now,  for  periodic  stresses, 
eiwt .  Equation  (14)  can  then  be 


both  a  and  S  have  time  dependence  of  the  form 
reduced  to 


Cm  = 


r-  2  m 

E  or  +  l  to - 

m  h 

m 


(15) 


From  the  wave  equation,  Eq.  (3),  it  was  shown  that  the  shear  modulus  operator 
(M)  could  be  reduced  to 

M  =  u  ♦  iuu'  .  (16) 

Therefore,  separating  real  and  imaginary  parts  of  Eq.  (15)  and  comparing  with 
Eq.  (16), 


so  that 


(17) 


(18) 


The  solution  given  by  Eq.  (9)  can  then  be  immediately  written  as 
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which  is  the  general  solution  for  the  Maxwell  model  of  a  visco-elastic  solid. 

To  evaluate  the  elastic  and  viscous  constants,  experiments  must  be  performed 
to  measure  the  acoustical  velocity  and  attenuation  at  a  given  frequency.  Equa¬ 
tion  (19),  with  the  assistance  of  Eq.  (17),  represents  two  simultaneous  equations 
with  two  unknowns  so  that  the  frequency-independent  parameters  Em  and  %  can  be 
derived. 

To  determine  the  applicability  of  this  model  at  high  and  low  frequencies, 

consider  the  following.  For  oj  »  E  /n 

m  m  ’ 
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Therefore 


k 


(23) 


which  is  the  solution  for  an  acoustic  wave  moving  in  an  elastic  medium.  There¬ 
fore,  for  this  model  the  medium  appears  to  be  elastic  for  velocity  considerations. 
The  attenuation  is  then  given  by 


T 


1 

C 


m  pc 

—  :  —  =  constant  . 

n  n 

m  m 


(24) 


For  the  low  frequency  case, 


u  << 


E 

m 


and 


but 


T  :  k 


pu)2 

TiT 


n  u> 

m 


m 


(25) 


(26) 


and 


V  *  -2^—  (27) 

m 

so 

T  ~  ~  .  (28) 


Thus  at  low  frequencies  both  the  attenuation  and  the  velocity  of  the  acoustic 
disturbance  have  frequency  dependence  of  the  form  »4). 
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Voigt  Model 

Figure  2  shows  the  mechanical  representation  of  a  visco-elastic  solid 
proposed  by  Voigt3.  In  this  configuration,  a  spring  with  elasticity  (Ev)  is 
in  parallel  with  a  dashpot  which  represents  the  viscosity  (nv)  •  As  in  the 
previous  model,  both  Ev  and  nv  are  assumed  to  be  frequency  independent. 


Figure  2.  Voigt’s  Mechanical 

Representation  of  Visco- 
Elastic  Solids. 


The  stress-strain  relationship  is  given  by 


o 

V 


(E  ♦  n 
v 


(29) 


Using  techniques  and  arguments  identical  to  those  for  the  Maxwell  model,  it  can 
be  shown  that  the  solutions  for  the  Voigt  model  are  given  by  the  equation 


and 


V  •  Ey 

u'  «  nv 


(30) 


(31) 
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cor  the  high  frequency  case,  i.e.,  where  Ev/ri is  small, 

T  =  k 


(32) 


which  is  the  same  form  as  the  low  frequency  limit  of  the  Maxwell  model,  for 
which  both  the  attenuation  and  the  velocity  vary  with  frequency  as 

For  the  low  frequency  case,  where  Ev/nvo)  is  large,  the  equations  approxi¬ 
mate  to 

n. 


T  :Tkr“ 

v 


(33) 


par 


so  that 


"“It 


(34) 


This  equation  states  that  the  velocity  remains  constant  while  the  attenuation 
varies  as  the  square  of  the  frequency. 


EXPANSION  TO  MORE  COMPLEX  MODELS 

Very  few  solids  behave  even  approximately  like  either  the  Maxwell  or  the 
Voigt  model.  However,  since  more  complicated  models  become  extremely  involved 
mathematically  and,  further,  since  models  specifying  more  parameters  require 
more  experimental  evidence  to  substantiate  or  "fit"  the  models,  using  a  single 
Maxwell  or  Voigt  element  is  often  a  convenient  method  of  obtaining  a  description 
of  a  visco-elastic  solid's  mechanical  properties  when  its  dynamic  behavior  is 
required  for  only  a  restricted  region  of  frequencies. 

For  most  materials,  the  use  of  single  elements,  although  mathematically 
more  desirable,  must  include  the  assumption  that  the  elastic  and  viscous 
coefficients  are  frequency  dependent.  To  eliminate  this  frequency  dependence, 
a  more  complex  model  becomes  necessary.  The  method  used  to  construct  a  more 
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complex  model  is  to  build  up  various  combinations  of  Voigt  and  Maxwell  ele¬ 
ments  in  series  and  parallel  until  the  frequency  dependence  of  the  constants 
is  removed. 

As  a  example,  consider  the  model  shown  in  Figure  3,  which  consists  of  a 
Maxwell  and  a  Voigt  model  in  series.  The  stress-strain  relationship  for  this 
configuration  must  be  reduced  to  that  described  by  Hook's  law,  which  is  given 
in  Eq.  (11).  From  earlier  arguments,  it  has  been  shown  that  the  Maxwell  and 
Voigt  elements  of  this  model  obey  the  equations 


> 


and 


_d_ 

dt 


E 

0  +  JE  0 
m  n  m 
in 


E  ± 
m  dt 


S 

m 


(14) 


o 


V 


(E, 


(29) 


respectively. 


Figure  S.  Mechanical  Representation  of 
the  4-Element  Model. 


To  combine  the  elements,  it  must  be  recognized  that  the  stress  across  the 
model  is  equally  divided  between  the  elements,  whereas  the  strain  is  the  sum  of 
the  extension  of  each  element.  Mathematically, 

o  =  ov  «  a 

s  =  sv  ♦  Sn  (35) 

so,  from  Eqs.  (10)  and  (29), 
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d  m 
■j -  a*  —  a 

dt  % 


=  E  j-  S 
m  dt  m 


(36a) 


o 


E  S  +  n  4r  S 
v  v  v  dt  v 


(36b) 


The  e1(i)t  dependence  cannot  be  substituted  directly  into  Eq.  (36)  because  each 
element  may  oscillate  at  frequencies  which  are  quite  different  from  the  driving 
frequency.  Thus,  before  the  time  dependence  can  be  considered,  the  generalized 
Hooke's  law  for  this  case  must  be  derived.  Rewriting  Eq.  (36b)  as 


E 

m 


o 


s  ♦  e  s 

v  m  dt  v 


and  adding  it  to  Eq.  (36a)  give 


(37) 


_d_ 

dt 


a  ♦ 


E  E 
m  v 


(S  -  S  )  ♦  E  S 
m  m  dt 


(38) 


Differentiating  Eq.  (38)  with  respect  to  time  gives 


E  E 
m  v 


_d_ 

dt 


(s  -  S  )  ♦  E 
m 


m 


S 


(39) 


Solving  Eq.  (36a)  for  (d/dtJS,,,,  substituting  into  Eq.  (39),  and  collecting 
terms  yield  the  relationship  between  the  time  derivatives  of  the  stress  and 
total  strain. 


dt 


a  * 


(E  E  E  \  ,  EE  ,2  EE 

d  a  +  JULa  x  E  d  S ♦JLSL 
\  nv  nv/dt  Vv  "dt2  nv 


dt 


S  . 


(40) 


Now,  for  periodic  driving  forces. 


„  _  „lUJt 

a  =  on  e. 

S  =  S°  elut 


(41) 


Substituting  Eq.  (41)  into  Eq.  (40)  and  simplifying  yield 
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But,  now 


M  =  u  +  iwM 


(42) 


so  that 


Rw  = 


0) 


U 


where 


(43) 


(44) 


To  investigate  the  behavior  of  this  model,  the  asymptotic  frequency  limits  must 
be  examined.  For  high  frequency, 
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but , 


so  that 


Lim 

t j )  -*large 

Lim 

a)  -►large 


M  =  Em 


Lim  (t) 
a)  -►large 


=  j  k  Roj 


Lim  (k)  =  £- 
<*»  -►large  m 


Lim  (t) 
a)  -►large 


1 

2 


=  constant  , 


(45) 


(46) 


(47) 


and,  for  low  frequency. 


Lim  (p)  =  const  x  m2 
co  -►small 


Lim  (Ruj) 
u)  -►small 


Lim  (t)  =  k 
w  ->small 


Lim  (k)  =  const  x  Vu  . 
oj  -►small 


(48) 


(49) 


Comparison  of  Eqs.  (46),  (47)  and  (19)  with  the  asymptotic  limit  solutions  of 
the  Maxwell  model  shows  that  both  models  have  the  same  form  of  frequency  de¬ 
pendence  in  the  limit  of  both  the  high  and  low  frequencies.  Therefore,  the 
Maxwell  element  predominates  the  frequency  dependence  of  this  model.  The 
limits,  of  course,  are  different  because  the  coefficients  are  highly  modulated 
by  the  Voigt  element. 
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EXTENSION  OF  THE  THEORY  TO  INCLUDE  ANISOTROPY 

Thus  far  in  the  theory,  only  the  solutions  for  an  isotropic  medium  have  been 
discussed.  It  is  known,  however,  that  for  many  materials,  including  naturally 
occurring  ice,  the  acoustic  properties  exhibit  a  dependence  on  the  direction  of 
the  propagation  of  the  wave  relative  to  the  crystal  axis,  i.e.,  the  media  are  not 
isotropic.  A  method  of  including  anisotropy  in  the  theory  follows. 


Stresses  and  strains  in  the  medium  can  be  specified  in  terms  of  the  resolved 
displacements  in  a  unit  cube  of  the  material.4  There  are  six  components  of  stress 
and  six  of  strain  which  can  be  written  in  terms  of  tensor  notation  as 

o  j j  and  j  , 

where  i,j  denote  the  x-,  y-,  z-  or  1-,  2-,  3-component  of  the  stress  or  strain. 

The  force  of  a  unit  cube  in  the  i^  direction  is  given  by 


3 o.  . 

=  -g-il  dxdydz  ,  (50) 

j 


where  repeated  indices  indicate  a  summation.  The  tensor  strains  are  defined  by 


S.  . 


i,j  =  1,2,3  , 


(51) 


where  the  uj 's  are  the 


displacements  of  the  body  along  the  coordinate  direction 


Utilizing  the  reduced  tensor  convention  for  the  stresses 


°1  '  °11;  °2  '  °22;  °3  ’  °33 

°4  *  °23  "  °32;  °5  '  °13  "  a31;  °6  *  °12  -  ”21 


(52) 


and,  similarly,  for  the  strains 

S1  •  SU:  S2  ■  S22;  S3  ■  S33 

21  "  S23  ’  S32*  r  -  SJ3  -  S3l:  /  “  S12  '  S21 


(53) 
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the  generalized  Hooke's  law  may  then  be  written 

°i  =  c  Sj  i,j  =  1,2. ...6  ,  (54) 

and  it  can  be  shown  that  for  the  elastic  constants 


cij  "  cji 


(55) 


there  are  a  maximum  of  21  elastic  constants  for  the  most  unsymmetrical  crystal. 
As  the  degree  of  symmetry  increases  the  number  of  elastic  constants  decreases 
until  there  are  only  two  for  an  isotropic  medium.  These  are  the  Lame'  constants 
X  and  p,  where 


X  +  2u  =  cu  =  C22  =  C33  , 


X  =  C.  „ 

=  c 

=  C„ 

12 

13 

23 

CJ 

II 

=  css 

C66 

and  all  other  constants  are 

zero. 

(56) 


To  include  anisotropy  in  the  theories  developed  in  this  report,  each  of  the 
elastic  constants  is  replaced  by  a  differential  operator  of  the  proper  form  for 
the  model  being  considered.  Therefore,  a  complete  solution  will  require  the 
determination  of  two  or  more  elastic  and  viscous  constants  for  each  operator 
used. 


SHORT  WAVELENGTH  SOLUTIONS 

GENERAL  DISCUSSION 

An  acoustic  wave  propagating  in  a  nonhomogeneous  medium  will  interact  with 
the  inhomogeneities  to  scatter  part  of  the  energy  out  of  the  sound  beam.  When 
the  wavelength  of  the  sound  is  large  compared  to  the  size  of  these  scattering 
centers,  th«;  interaction  is  minimal  and  the  energy  loss  is  small  compared  to 
that  caused  by  classical  dissipation.  For  this  situation,  the  acoustic  attenu¬ 
ation  can  be  predicted  by  the  long  wavelength  theories  discussed  in  the  previous 
section.  However,  when  the  size  of  the  scattering  center  is  a  significant 
portion  of  a  wavelength,  the  contribution  to  the  attenuation  due  to  scattering 
processes  cannot  be  neglected. 

To  determine  the  attenuation  due  to  scattering,  consider  an  isotropic 
solid  elastic  medium  containing  a  concentration  of  scatterers.  The  inhomo¬ 
geneities  in  the  medium  can  be  mathematically  characterized  by  a  generalized 
function, 
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D(S, I ,a,r . . .)  =  S(si)I(si,sj)G(a,r)...  ,  (57) 

where  S(sj)  is  a  structure  distribution  function  which  describes  the  relative 
importance  of  the  various  types  of  scatterers  in  the  medium,  I(s^,sp  is  a 
distribution  function  which  describes  the  interaction  between  scatterers  and 
multiple  scattering  of  the  sound,  and  G(a,?)  describes  the  distribution  of 
the  size  and  position,  etc.  of  the  scatterers. 


Assuming  a  plane  wave  moving  in  the  positive  x-direction,  the  intensity  of 
the  sound  is  given  by 


I 


-2ax 


(58) 


where  a  is  the  attenuation  constant,  x  is  the  penetration  distance  into  the 
material  and  I0  is  the  intensity  of  the  original  sound  wave.  If  the  energy 
losses  are  all  considered  to  be  due  to  scattering  processes,  Eq.  (58)  may 
be  rewritten 


I  =  I0  e'Fx  ,  (59) 

where  T  is  the  generalized  scattering  cross  section  and  is  governed  by  the  size, 
shape,  position,  etc.  of  the  scatterers  in  the  medium.  In  terms  of  Eq.  (57),  T 
may  be  written 


r  *£//••  .D(S,I,G,r, . . .)  YiUJdrda... 

i.j 


(60) 


r  a 


where  nowy^(a)  is  the  scattering  cross  section  of  a  single  independent  scatterer. 
The  attenuation  is  then  given  by 


.D(S,I,G,r,...)  Yi(a)drda...  .  (61) 

lj  -► 

J  r  a 

This  equation  states  that  the  attenuation  is  determined  by  the  scattering 
cross  section  and  the  distribution  of  the  scatterers.  The  cross  section  is  de¬ 
pendent  upon  the  interaction  of  the  wave  with  each  individual  scatterer  while 
the  distribution  function  is  the  probability  of  the  wave  encountering  that 
scatterer  in  the  medium. 


SCATTERING  OF  A  PLANE  WAVE  FROM  A  SPHERICAL  OBSTACLE 

Consider  the  problem  shown  in  Figure  4.  An  acoustic  plane  wave  propagating 
in  an  isotropic  elastic  solid  encounters  a  spherical  scatterer  which,  for  present 
discussion,  is  also  assumed  to  be  elastic.  Both  media  are  characterized  by  their 
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densities  (p^),  shear  moduli  (Uj)  and  Lame'  constants  C^i) -  Elasticity  theory 
defines  the  compressicnal  and  transverse  velocities  as 


/ 


Xi  +  2yi 


pi 


(62) 


respectively,  where  u>  is  the  circular  frequency  and  ki  and  Kj  are  the  wave- 
numbers.  The  displacement  (u)  can  be  written  in  terms  of  the  displacement 
potentials  0  and  ip  as 

u  =  grad  <p  *  curfl.  $  .  (63) 


INCIDENT 

WAVE 


SCATTERER 


SCATTERED 

WAVE 


Figure  4.  Representation  of  the  Problem. 

In  spherical  coordinates, 

u  =  -V$  +  VxVx(ri|/)  ,  (64) 

which  is  valid  for  spherically  symmetric  waves.  The  displacement  potentials  in 
spherical  coordinates  are  of  the  form 

00 

E  Cnfjn(Pr)  PB(cos0)  eiwt  ,  (65) 

m-o 
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and  must  obey  the  wave  equations 


where 


V2<J> 


i  a2<fr 

CL*  It* 


V2\J> 


1  32\|/ 
C?2  W 


(66) 


fm(pr)  -  the  spherical  Bessel  function  of  order  m 

p  :=  either  k  or  K,  depending  upon  the  wave  considered 

Pm  (cos0)  =  the  Legendre  polynomial  of  order  m. 

'/.lien  the  incident  plane  wave  strikes  the  scattering  center,  some  of  the  energy 
is  scattered  away  and  some  is  transmitted  through  the  scatterer.  Since  all  the 
waves  must  obey  the  same  wave  equations,  the  wave  potentials,  neglecting  the 
time  dependence,  ray  be  written  as 

00 

<j>  Incident  *  £  I  i  (k.r)  P  (cos8) 

cm  mi  m 

m*o 

00 

1||  Incident  =  £  I  j  (K,r)  P  (cos0)  (67) 

sm  m  i  m 

m«o 


00 

<J)  Scattered  *  £  A  h  (k.r)  P  (cos0) 

m  m  i  m 
m*o 


ip  Scattered  5  J]  B  h  (K.r)  P  (cos0)  ,  (68) 

^  m  m  1  m 
m«o 


<p  Transmitted  -  £  c.  L  (M)  P„  (cos0) 
m»o  ■  1  1  m 

00 

Transmitted  ■  £  j  (l^r)  P  (cos0)  (69) 

m«o 
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where  r  and  9  are  spherical  coordinates.  The  spherical  Bessel  functions  of  the 
first  kind,  jm(kr),  and  third  kind,  hm(kr),  are  used  to  assure  that  the  scattered 
wave  is  propagating  away  from  the  scatterer  located  at  the  origin  of  the  coordi¬ 
nate  system.  In  order  for  the  incident  wave  to  be  plane  and  of  unit  amplitude, 
Icm  and  Ism  must  be  of  the  form5 


/  .  V  X 

I  =  (2m  +  1)  ,  (70) 

m  p 

where  p  is  k  and  K,  respectively.  The  coefficients  Am,  Bm,  Cm  and  Dm  are  deter¬ 
mined  by  the  boundary  conditions  at  the  scattering  surface.  For  a  solid-solid 
interface  with  no  slippage  these  boundary  conditions  are 

(a)  Radial  displacement  continuous 

(b)  Tangential  displacement  continuous 

(c)  Radial  stress  continuous 

(d)  Tangential  stress  continuous. 

(When  the  scattering  spherj  is  fluid  filled,  slippage  is  allowed  2nd  the  require¬ 
ment  that  the  tangential  displacement  be  continuous  is  no  longer  valid.)  At  the 
boundary,  r  =  a,  these  conditions  can  be  written  mathematically  as 


(a) 

ur(inc) 

+  ur(scatt) 

ur(trans) 

(b) 

u0(inc) 

+  u0(scatt) 

u0(trans) 

(c) 

°rr(inc) 

+  °rr(scatt) 

°rr(trans) 

(d) 

°9r(inc) 

+  ar0(scatt) 

°r0(trans) 

In  terms  of  the  displacement  potentials,  the  stress  and  displacement  components 
of  Eq.  (71)  are6 
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where 


n  = 


1  a 


sin  0  3^ 


(sin9le) 


(73) 


The  scattering  cross  section  is  determined  by  calculating  the  ratio  of  the 
time  rate  at  which  energy  is  scattered  out  of  the  wave  by  an  obstacle  of  radius 
"a"  to  the  total  energy  of  the  incident  wave.  The  scattered  energy  is  equal  to 
the  energy  being  carried  away  by  the  scattered  wave  across  a  spherical  surface 
of  radius  "b"  >  "a".  Love4  gives  this  scattered  energy  as 


'scattered 


If 

A 


3u 


3u 


Jxr  5T  +  °yr  3t 


+  o 


3u 
_ 2 

ZT  3t 


dA 


scattered 


(74) 


where  axr,  a  and  ozr  are  the  stress  components  acting  in  the  three  rectangular 
axes  on  a  surface  normal  to  the  radius  vector  r,  and  the  u's  are  displacements. 
Both  the  a's  and  the  u's  are  generally  complex,  so  car*  must  be  exercised  to 
assure  that  the  final  expression  for  the  scattered  energy  is  real.  Assuming  a 
time  dependence  of  eiwt  for  the  displacement  potentials  and  using  the  spherical 
symmetry  of  the  scattered  wave,  Eq.  (74)  can  be  rewritten  in  spherical  coordinates 
as 


E 


scattered 


—  /  1 1°  u  +0o  uQ  +  o  u  j 

J  ||  rr  r  0r  0  zr  z  J 

r*b 


•  °rrur  +O0ru0  +  a 


""•II 


2rb2  sin0d0 

scattered 

wave 


(75) 


where  all  the  terms  are  defined  by  Eq.  (72)  and  the  asterisk  denotes  the  complex 
conjugate.  The  integral  in  this  equation  can  be  evaluated  by  substituting  Eqs. 
(68)  and  (72)  and  using  the  identity  which  is  valid  for  any  two  arbitrary 
functions  f  and  g. 


fjfr f)  g  sin0d0  •/JfW  sin0d0  -  -fj  ||||  sin0d0  ,  (76) 

to  determine  the  recursion  relations  and  the  orthogonality  of  the  Bessel  functions 
and  the  Legendre  polynomials.  Then,  the  energy  lost  from  the  incident  wave  due 
to  scattering  from  the  spherical  scatterer  is 
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Escattered  =  ^irp^ai 


OO 


E 

m=n 


l 


2m  +  1 


r-  |A  I 

kj  1 


+ 


m(m+l) 


(77) 


where  Am  and  B^,  are  the  magnitudes  of  the  scattered  wave  potentials  and  must  be 
evaluated  at  the  scattering  surface  by  the  boundary  conditions. 


The  total  energy  of  the  wave  is  determined  by  calculating  the  energy  in  the 
incident  plane  wave  that  is  being  transported  through  a  unit  area  normal  to  the 
propagatioi  direction.  The  displacement  of  a  plane  compressional  wave  of  unit 
amplitude  traveling  in  the  positive  z-direction  can  be  written  as 


u  =  z 


ei(o)t  -  kjZ) 


(78) 


From  elasticity  theory,  the  stress  can  be  calculated  from  the  relation 


where 


a..  «  X(V«u)  6 

ij  13 


1  i=j 

0  i/j 


/3u.  3u.\ 


(79) 


For  the  wave  described,  =  0  for  all  ij  except 


a  =  .  i  J.  ei(u)t-klz) 
zz  ki 


The  energy  flux  through  a  unit  area  is  then  given  by 


p.u)J 


(80) 


*  /  .  KI ^ 

E.  *  (o  u*  ♦  0  u*  ♦  a  u*)  -  (0  *u  +  a  *u  +  0  *u  )  >  =  -r— 
inc  2  \v  xz  x  yz  y  zz  z;  v  xz  x  yz  y  zz  z J  j  kj 

A  similar  calculation  yields  the  energy  flux  for  a  plane  shear  wave  to  be 

,s 


(81) 


"inc 


(82) 


The  scattering  cross  section  is  defined  as 
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Y  = 


'scattered 


'total 


(83) 


The  scattering  cross  section  for  a  compressional  plane  wave  encountering  a 
spherical  obstacle  is  then 

Y shear  =  4ir  E  2m  +  1  I  ^  AnJ  *  +  M  i  l  (84) 

m=o  [  1  J 

and  that  for  a  transverse  plane  wave  is 

Yco«p- 4,r  f  jsrrfr-  Ib.1’1  •  (*=) 

m=o  (.  1  J 

If  only  spherical  scatterers  are  present  in  the  medium,  Eqs.  (84)  and  (85) 
can  be  substituted  into  Eq.  (61)  to  determine  the  compressional  and  shear  attenu¬ 
ation.  Of  course,  before  this  can  be  accomplished  the  distribution  function  for 
the  material  must  be  known.  Also,  if  other  scatterer  configurations  are  present, 
scattering  cross  sections  for  each  type  must  be  calculated  using  techniques 
similar  to  those  above. 


APPLICATION  OF  THE  THEORIES  TO  SEA  ICE 

SCATTERING  ATTENUATION  APPROXIMATION 

When  sea  ice  freezes,  the  brine  in  the  water  concentrates  in  small  pockets. 
The  result  is  a  system  composed  of  frozen  water  with  imbedded  scattering  centers 
filled  with  a  brine  fluid  which  is  highly  saline.  Although  these  inclusions  can 
have  various  -geometries,  the  vast  majority  of  those  observed  in  annual  sea  ice 
by  the  authors  of  this  report  were  spheroidal.  This  same  situation  may  not  nec¬ 
essarily  hold  for  other  types  of  sea  ice,  e.g.,  multi-year  ice,  but  the  acoustic 
and  ‘petrographic  data  necessary  for  evaluation  of  other  types  of  sea  ice  are 
presently  unavailable.  Therefore,  for  the  purposes  of  this  calculation,  it  will 
be  assumed  that  all  of  the  scattering  centers  are  brine-filled  spheres. 

Consider  the  application  of  the  previous  theory  to  the  case  of  a  plane  wave 
being  scattered  from  a  liquid-filled  sphere.  The  amplitudes  of  the  scattered 
and  transmitted  waves  are  evaluated  by  solving  the  boundary  conditions  at  the 
scattering  surface  for  a  wave  striking  a  liquid-solid  interface.  By  inserting 
Eqs.  (67)  through  (69)  into  Eqs.  (71a),  (71c)  and  (71d)  and  evaluating  at  r  =  a, 
three  simultaneous  equations  result.  They  are 
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and 


■cWA')-1.'  j..l(kla)-Is.['(m*1)  W)] 

'  -*.l,h»(kla)  -klah^l(kla))  *  “.I"'"*11  VKla)1 

*  Cni"Vk2a)  -  k2a  Vl(k2a»  • 


cm 


m(l-m)  ♦ 


h*** 


jBtki»)  -  2kia  j^dcja) 


♦  Isii[m(».l)]  [(n-1)  j^CKjSj-Kja  j^jO^a)] 

K  2a2l 

■(1-m)  — I  h^a)  -  2kja  h^O^a) 


=  -A 


Bm[«Cm*l)]  [(■-!)  hBCK1a)  -  l^a  h^O^a)] 


P2  S2a2 


+  p,  2  Cm  jm(k2a)  * 


(86a) 


(86b) 


W*-1)  j.(kia)  -kia  Wk!a)J 

♦  ‘J-^)  -  v  vi<kia>] 

*  -AbI(«-1)  hB(kxa)  -kja  h^jCkja)] 

•  B.  [(1_“2  +  ‘1r")h«(K1*)  _  Kla  WKla)]  •  (86c) 


These  three  equations  must  be  solved  for  the  amplitudes  Ag  and  Bg,  to  evaluate  the 
scattering  cross  section. 

In  Eq.  (86),  1^,  and  Ism  are  the  amplitudes  of  the  compress ional  wave  and 
the  shear  incident  wave,  respectively.  To  consider  the  propagation  of  a  shear 
wave  in  the  medium,  the  amplitude  of  the  compressional  wave  must  vanish,  i.e., 

Icm  *  0.  It  was  shown  earlier  in  Eq.  (70)  that  Isa  is  given  by 

rsm "  iq- (“i)*^1(2^1)  *  (87) 
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Using  the  values  for  the  velocities  and  densities  of  annual  sea  ice  from  the  ex¬ 
periment  described  in  Appendix  A  and  the  computer  program  SCTTR  given  in  Appendix 
B,  we  calculated  the  product  of  the  scattering  cross  section  and  the  frequency 
squared  (yv2)  vs  the  circumference-to-wavelength  ratio  of  the  scattering  sphere. 
The  results  of  this  computation  are  shown  in  Figure  5.  In  the  frequency  region 
for  which  ka«l,  i.e.,  for  wavelengths  much  larger  than  the  scatterer  radius, 
the  scattering  cross  section  is  of  the  form 

’'shear  '  “6  ■  (88) 

which  is  the  result  expected  from  standard  Rayleigh7  scattering.  For  frequencies 
such  that  ka»l,  the  scattering  cross  section  becomes  of  the  form 


Y  ~  \T’°  aJ* 

'shear  *  (89) 

ka>10 

In  general,  the  shear  wave  scattering  cross  section  for  ka<l  and  ka>15  can  be 
approximated  by 


v4a6 

rshcar  *  4.151  x  10'9  ♦  2.35  x  10“5v2,2a2,2 


dB/cm 


(90) 


where  the  frequency  (v)  is  in  MHz  and  the  scatterer  radius  (a)  is  in  cm.  Thesi 
expressions  do  not  describe  the  erratic  behavior  of  the  scattering  cross  section 
where  ka~l.  The  fluctuations  in  this  region  are  due  to  bubble  resonance  and 
affect  the  cross  section  over  a  small  ka  span.  However,  for  efficient  propaga¬ 
tion  of  the  wave  through  the  medium,  frequencies  with  wavelengths  approximately 
the  same  size  as  the  scattering  centers  should  be  avoided. 

For  the  case  of  a  compressional  wave  incident  on  the  bubble,  Ism  *  0  and 

^  (-i)*+1(2»+l)  .  (91) 

The  solution  for  the  scattering  cross  section  is  determined  following  the  same 
procedure  described  for  the  transverse  wave.  The  results  of  this  computation 
are  also  plotted  in  Figure  5  as  a  function  of  ka.  For  values  of  ka<l  and  ka>15, 
the  scattering  cross  section  can  be  approximated  by  the  equation 


4  6 

_ v  a _ _ _ _ 

Yco«p  2.163  x  10-7  ♦  0.017  v4  a* 


dB/cm 


(92) 
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Figure  5.  Results  of  Scattering  Cross  Section  Calculation. 
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where  v  is  in  MHz  and  a  is  in  cm.  As  in  the  solution  to  the  shear  wave  case, 
at  low  values  of  ka  the  compressional  wave  scattering  cross  section  is  dominated 
by  Rayleigh  scattering.  For  large  values  of  ka,  however,  the  cross  section  be¬ 
comes  a  constant  for  the  compressional  wave  while  the  transverse  wave  approaches 
infinity.  This  result  indicates  that  when  the  wavelength  becomes  short  the  rel¬ 
ative  curvature  of  the  scatterer  becomes  very  large,  so  that  the  scattering  prob¬ 
lem  approaches  that  of  the  transmission  and  reflection  of  an  acoustic  wave  from 
a  liquid  layer  in  the  ice  medium.  In  this  case,  changes  in  frequency  no  longer 
affect  the  relative  transmission  of  the  compressional  wave  energy  into  the  fluid, 
and  the  reflectivity,  i.e.,  the  scattered  energy,  remains  constant.  Conversely, 
if  a  shear  wave  strikes  this  boundary,  it  is  totally  reflected,  i.e.,  the  scat¬ 
tering  cross  section  is  infinite,  because  the  fluid  medium  will  not  support 
shear  waves. 

To  evaluate  the  attenuation  due  to  scattering,  the  distribution  functions  in 
Eq.  (61)  must  be  either  known  or  assumed.  Because  the  necessary  data  are  presently 
unavailable,  and  to  minimize  the  complexity  of  the  calculation,  the  following 
assumptions  have  been  made. 

(a)  No  interaction  exists  between  the  scattering  centers,  and  the  wave, 
once  scattered,  will  not  be  re-scattered  back  into  the  sound  beam. 

Since  second  order  scattering  should  have  only  a  small  effect,  par¬ 
ticularly  in  the  Rayleigh  scattering  region,  this  assumption  should 
not  greatly  modify  the  calculation  accuracy. 

(b)  There  are  N  scattering  centers  per  unit  volume  and  the  centers  are 
uniformly  distributed  throughout  the  medium.  From  observations 

in  the  laboratory,  this  assumption  appears  to  be  relatively  valid. 

fc)  All  of  the  scattering  centers  are  of  the  same  size,  with  a  radius 
aQ.  This  is  probably  an  invalid  assumption  because  it  is  known 
that  a  wide  variety  of  brine  pocket  sizes  is  present  in  the 
medium.  However,  if  radius  aQ  is  selected  as  the  average  bub¬ 
ble  size,  then  it  can  be  argued  that  only  the  larger  bubbles  in 
the  distribution  will  significantly  affect  the  scattering.  This 
assumption  will  make  the  calculated  attenuation  less  than  the 
actual  attenuation  that  would  be  expected  to  occur  during  an 
experiment . 

Using  these  assumptions,  the  attenuation  due  to  scattering  can  be  reduced  to 

||| 

Scattering  "  7  ^ao^  *  (93) 

which  predicts  the  minimum  attenuation  to  be  expected  from  this  mechanism. 

Laboratory  measurements  of  void  size  and  concentration  yield  radii  on  the 
order  of  0.065  cm  and  concentrations  on  the  order  of  25  holes  per  cm3.  Using 
the  curves  of  Rigure  5  and  the  relation 
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(94) 


where  v  is  the  frequency  in  Hz,  c  is  the  velocity  in  cm/sec,  and  a  is  the  radius 
in  cm,  the  attenuation  due  to  scattering  at  500  kHz  is  calculated  as 


“scattering  comp  =  25  dB/meter 

and 

“scattering  shear  =  350  dB/meter  .  (95) 

At  100  kHz 

“scattering  comp  =  0*045  dB/meter 
and 

“scattering  shear  =  2-25  dB/meter  .  (96) 

In  general,  for  ka<l  and  ka>15 


v 

“shear  ~  -5  -4  2  2 

Snear  4.403  x  10  ♦  6.097  x  10  V*** 


dB/m 


and  (97) 

4 

a  -  - - — - s — t  dB/m  , 

COrap  2.294  x  10"3  ♦  3.253  x  10'3  v4 


where  v  is  in  MHz.  For  values  of  l<ka<15,  the  approximations  of  Eq.  (92)  are 
still  more  or  less  valid  for  order  of  magnitude  calculations,  but  will 
invariably  yield  predictions  which  are  significantly  low. 


TOTAL  ATTENUATION  APPROXIMATION 

As  was  discussed  in  the  introduction,  the  attenuation  in  sea  ice  is 
composed  of  two  parts, 

a  =  “classical  *  “scattering  *  (98) 

Before  the  total  attenuation  can  be  estimated,  the  classical  attenuation  must  be 
considered.  The  experiment  outlined  in  Appendix  A  gives  the  results  necessary 
for  an  approximation  to  the  solution  of  thu  two-element  classical  models.  During 
this  experiment,  the  total  attenuation  and  velocity  were  measured  in  situ  15  cm 
below  the  surface  in  annual  sea  ice.  At  100  kHz,  the  results  of  tKFse  measure¬ 
ments  gave  a  longitudinal  wave  velocity  of  3790  m/sec,  a  shear  wave  velocity  of 
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1923  m/sec,  a  density  of  0.937  g/cm3,  and  a  generalized  attenuation  of  0.18 
dB/cm.  Using  these  values,  and  assuming  the  elastic  and  viscous  constants  to 
be  frequency  independent,  we  solved  both  the  Voigt  and  the  Maxwell  model  as  a 
function  of  frequency.  The  results  are  shown  in  Figure  6,  along  with  the  scat¬ 
tering  attenuation.  From  this  representation  it  can  be  concluded  that,  for 
frequencies  above  approximately  5  kHz,  the  Maxwell  model  fails  because,  from 
the  physics  of  the  problem,  there  is  no  frequency  region  over  which  the  attenu¬ 
ation  remains  constant.  For  frequencies  lower  than  1  MHz,  these  curves  show 
that  both  the  attenuation  predicted  by  the  Voigt  model  and  that  predicted  by 
the  scattering  calculations  can  be  found  by  utilizing  their  low  frequency  ap¬ 
proximations.  Thus,  if  it  is  assumed  that  the  total  attenuation  is  the  sum  of 
the  results  of  the  Voigt  model  and  the  scattering  calculations,  then  the  attenu¬ 
ation  will  increase  as  v2  for  frequencies  below  approximately  250  kHz  and  as  v1* 
for  higher  frequencies  because  of  the  dominance  of  the  scattering  contribution. 
The  total  attenuation  curve  for  the  shear  wave  is  shown  in  Figure  7. 


DEPTH  DEPENDENCE  OF  SHEAR  WAVE  PROPAGATION 

Recent  experimental  data8  have  demonstrated  that  for  annual  sea  ice  the 
velocity  of  longitudinal  waves  exhibits  dependence  on  the  depth  into  the  medi¬ 
um  because  of  variations  in  density,  salinity  and  temperature.  It  is  reasonable 
to  assume  that  these  same  parameters  will  similarly  affect  the  propagation  of 
the  shear  wave.  Therefore,  it  is  desirable  to  prognosticate  the  form  of  the 
depth  dependence.  To  accomplish  this  calculation,  it  is  assumed  that  the  mea¬ 
surement  frequency  is  sufficient  to  ignore  the  viscous  terms  in  the  propagation 
equation,  or,  equivalently,  that  the  viscous  terms  do  not  greatly  affect  the 
velocity.  Thus,  a  near-elastic  medium  is  assumed,  for  which,  according  to 
elasticity  theory,  the  longitudinal  and  transverse  wave  velocities  are  given  by 


r  /E  (1  -  e)  l\*2u 

CL  “V  P  Cl  ♦  e)  (1  -  2e)  *V  P 

(99) 

•t  'vfirSy  'VF  • 

where 

Cl  is  the  velocity  of  the  longitudinal  wave 
Cj  is  the  velocity  of  the  shear  wave 
e  is  Poisson's  ratio 
E  is  Young's  modulus 
U  is  the  shear  modulus 
X  is  the  Lame"  constant 
p  is  the  density. 
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Figure  7.  Predicted  Attenuation  in  Annual  Sea  Ice. 

Langleben*  has  empirically  related  Yeung's  modulus  of  sea  ice  to  the  brine 
content  (v)  by  the  equation 

E  =  (10.0  -  35.1  v)  x  1010  dyne/ cm2  .  (100) 

The  brine  content  in  this  equation  can  be  determined  from  the  salinity  (S)  and 
the  temperature  (0)  by  the  relations10 

V  «  S  (-§-  ~  2.28 j  -  0.5°  <  0  <  -  2.6° 

v  -  s  (45q917+0.93o)  -  2.6*  5  0  i-8.2#  ,  (101) 

v  »  S  (—q— 5 ♦  1.189)  -  8.2°  <  0  <- 27.9° 
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where 


S  is  the  salinity  in  parts  per  thousand 
v  is  the  brine  content  in  parts  per  thousand 
0  is  the  absolute  value  of  the  temperature  in  °C. 

Using  these  equations  and  the  data  in  Ref.  8,  the  depth  dependence  of  the 
shear  wave  velocity  has  been  calculated  and  is  shown  in  Figure  8.  The  calcula¬ 
tion  shows  that  the  shear  wave  velocity  decreases  continuously  with  depth  and 
may  approach  zero  at  the  ice-water  interface  (a  depth  of  14S  cm  for  this  experi¬ 
ment).  This  decrease  in  velocity  can  be  related  to  the  vertical  distribution  of 
the  liquid  brine  present  in  the  ice  cover.  The  result,  then,  indicates  that  the 
boundary  might  not  provide  a  good  reflecting  surface  for  the  transverse  acoustic 
mode. 


Unfortunately,  sufficient  data  do  not  exist  to  similarly  predict  the  depth 
dependence  for  attenuation.  For  this  determination,  the  form  of  both  the  elastic 
and  viscous  coefficients  must  be  known.  However,  if  a  near-elastic  medium  is 
assumed,  it  can  be  shown  that  for  a  Voigt  solid 


T 


pnoj2 

2C* 


(102) 


Therefore,  provided  the  depth  dependence  of  the  viscosity  (n)  does  not  decrease 
faster  than  C3,  the  attenuation  will  become  infinite  as  the  shear  velocity  goes 
to  zero  at  the  boundary. 


CONCLUSIONS 

The  attenuation  of  an  acoustic  wave  in  any  material  is  a  result  of  the 
combination  of  classical  dissipation  and  scattering  effects.  For  sea  ice  it 
was  shown  that  the  Maxwell  model  for  classical  attenuation  does  not  apply  to 
the  dynamic  case  because  it  predicts  constant  attenuation  at  frequencies 
above  a  few  kHz.  This  conclusion  probably  can  be  extended  to  more  complex 
models  that  have  Maxwell  elements  in  series  because  it  was  shown  that  when¬ 
ever  this  component  is  in  that  position  it  dominates  the  model.  Therefore, 
if  more  complex  models  are  required  to  predict  experimental  values  using  fre¬ 
quency  independent  parameters,  combinations  of  Voigt  and  Maxwell  elements  in 
parallel  will  be  required. 

A  minimum  attenuation  due  to  scattering  was  derived  by  solving  for  the 
scattering  cross  section  of  a  single  liquid-filled  spherical  scatterer, 
assuming  a  uniform  bubble  population  and  size.  It  was  shown  that  for  the 
frequencies  of  interest,  Rayleigh  scattering  predominates  where  the  atten¬ 
uation  is  predicted  to  be  proportional  to  the  fourth  power  of  the  fre¬ 
quency.  Assuming  a  dissipation  that  could  be  predicted  by  a  single  Voigt 
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Figure  8 
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element,  we  derived  a  frequency-dependent  curve  composed  of  the  attenuation 
predicted  from  both  classical  and  scattering  sources.  This  calculation  was 
based  on  measured  values  of  sea  ice  velocity  and  attenuation  at  100  kHz. 

The  results  indicate  that  for  frequencies  above  250  kHz  scattering  dominates, 
and  the  attenuation  for  the  shear  wave  will  be  in  excess  of  100  dB/meter. 

Therefore,  based  upon  these  predictions,  all  experiments  for  which  an 
acoustic  shear  wave  is  to  be  used  and  for  which  penetration  into  the  medium 
is  desired  (such  as  sea-ice  thickness  measurements)  should  use  frequencies 
of  20  kHz  or  lower  depending  upon  the  distance  the  sound  wave  is  expected 
to  travel.  For  velocities  on  the  order  of  2000  m/sec,  these  frequencies  will 
have  wavelengths  of  10  cm  or  longer.  At  these  wavelengths,  resolution  is 
very  poor. 

The  probability  that  attenuation  is  a  continuously  varying  parameter 
depending  on  the  distance  from  the  ice  surface  was  discussed.  Insufficient 
experimental  data  presently  exist  to  verify  or  deny  this  assumption.  How¬ 
ever,  since  the  temperature  of  sea  ice  varies  with  the  depth  from  the  surface, 
higher  concentrations  of  liquid  brine  could  result  which  would  contribute  to 
a  decrease  in  shear  wave  velocity  and  an  increase  in  the  expected  attenuation. 
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APPENDIX  A 

DESCRIPTION  OF  IN  SITU  SHEAR  WAVE  EXPERIMENT 

During  March  1973,  an  extensive  experiment  was  performed  in  the  Chukchi 
Sea  near  Pt.  Barrow,  Alaska,  to  determine  the  acoustic  propagation  character¬ 
istics  of  annual  sea  ice.  Although  these  experiments  succeeded  in  measuring 
compressional  velocity  profiles  for  the  first  time,  similar  profiles  for  the 
shear  wave  were  not  obtained.  The  principal  difficulty  in  making  shear  wave 
measurements  is  obtaining  a  reliable  bond  between  the  medium  and  the  shear  wave 
generating  transducers.  This  difficulty  is  largely  due  to  the  "leaching"  of 
the  brine,  which  generates  a  highly  saline  liquid  surface  layer. 

Washing  the  bonding  surface  with  fresh  water  has  been  partially  success¬ 
ful  in  both  laboratory  and  in  situ  tests.  The  experiment  shown  in  Figure  A1 
was  attempted  during  the  field  tests  discussed  above.  Channels  were  cut  in 
the  ice  canopy  to  a  depth  of  ~45  cm,  and  located  a  known  distance  apart. 

The  sides  of  the  channel  were  washed  with  fresh  water  to  eliminate  the  saline 
boundary  layer  and  shear-sensitive  transducers  were  quickly  coupled  to  the 
ice  at  a  depth  of  15  cm  on  diametrically  opposite  sides  of  the  test  section. 

The  electronics  consisted  of  a  pulse  timing  generator  which  furnished  PRF 
and  pulse  length  information  to  a  pulsed  power  oscillator  and  a  trigger  sig¬ 
nal  to  the  oscilloscope.  The  oscillator,  in  turn,  provided  a  cw  pulse  of 
the  desired  frequency,  pulse  length  and  pulse  rate  to  the  transmitting  trans¬ 
ducer  coupled  to  the  ice.  The  transmitted  acoustic  signal  was  received  by  the 
receiving  transducer  coupled  to  the  opposite  side  of  the  test  section,  ampli¬ 
fied  and  transmitted  to  the  oscilloscope.  The  time  of  flight  of  the  acoustic 
pulse  in  the  test  section  was  determined  by  measuring  the  oscilloscope  sweep 
time  from  the  origin  to  the  first  received  signal  by  means  of  a  calibrated 
time  delay  and  then  correcting  for  electronic  time  losses.  Thus,  knowing  the 
specimen  length,  the  velocity  could  be  calculated.  This  method  is  accurate  to 
approximately  1-1/2%,  depending  on  how  precisely  the  length  of  the  sample  is 
known. 

This  technique  was  used  in  a  series  of  measurements  during  which  the  test 
section  was  incrementally  decreased  and  the  process  repeated.  These  measure¬ 
ments  yielded  an  s  erage  attenuation  of  18  dB/meter  and  an  average  velocity  of 
1923, m/sec  at  100  kHz.  The  velocity  is  shown  as  an  experimental  point  in 
Figure  8.  Further  laboratory  measurements  of  shear  wave  velocity  on  sea  ice 
samples  using  techniques  similar  to  those  of  the  ui  situ  experiment  showed  an 
average  velocity  of  1870  m/sec,  which  agrees  very  well  with  the  average 
predicted  profile  velocity  of  1874  m/sec. 
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Figure  Al.  Diagram  of  the  in  Situ  Experiment. 
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APPENDIX  B 

COMPUTER  PROGRAMS 


PROGRAM  SSCTRt INPUT .OUTPUT, TAPE5»IAPUT,TAPE6«CUTPUT) 

THIS  PROGRAM  CALCULATES  THE  SCATTERING  CROSS  SECTION  TIMES  THE 
FREQUENCY  SQUAREO.  THE  VALUE  IS  CALCUcATEO  FOR  A  SHEAR  ®LANE 
NAVE  IN  A  SOLIO  STRIKING  A  FLUIO  FILLED  SPHERE.  THE  CASE  OF  A 
COMPRESS IONAL  NAVE  STRIKING  THE  SPHERE  IS  ALSO  CALCULATED.  THE 
VALUE  IS  CALCULATEO  AS  A  FUNCTION  OF  THE  OIMENSIONLESS  PARANETFR 
KA  MHERE  K*2*PI*FREQ/ VELOCITY  AND  A  IS  THE  SPHERE  RADIUS. 

THE  CALCULATION  IS  NAOE  IN  THE  C .G.S. SYSTEM.  THE  SCATTERING 
CROSS  SECTION  HAS  UNITS  OF  09-CN**2.  CS1  IS  THE  SHEAR  VELOCITY 
ANO  CL1  IS  THE  CONPRESSIONAL  VELOCITY  IN  THE  SOLIO.  CL2  IS  THE 
COMPRESS IONAL  VELOCITY  IN  THE  SCATTE»ER.  RHO  IS  THE  RATIO  OF  THE 
OENSITY  OF  THE  SPHERE  TO  THE  OENSITY  OF  THE  SOLIO. 

READ(5,*9>  CSl.CLl.CL2.RH0 
*9  FORMAT (3F8.0.F6.3) 

NRITE (6,51) 

51  FORMAT U2H  <A  GAMMA, SHEAR  GAMMA, COHPI 

1  READC5.50)  AKA 
51  FORMAT (FT. 21 
IFCE0F.5I  7,8 
8  CONTINUE 

PI-3.1A159265A 

FIRST  THE  SHEAR  SCATTERING  IS  CALCULATED. 

9KS1*AKA 

BKL1*BKS1#CS1/CL1 
BKL2*9KS1*CS1/CL2 
SIM»CS1/«2.*PII 

IF  CIM*0.  THEN  A  SHEAR  NAVE  IS  INCIOENT. 

CIMM. 

CALL  COEF (BKS1,9KL1,9KL2, SIM, C IN, RHO, GAMAS) 

GAMAS*GAMAS*8KS1 

GAMAS  IS  THE  SHEAR  SCATTERING  CROSS  SECTI0N*FREQ**2. 

NON  THE  CONPRESSIONAL  SCATTERING  IS  CALCULATEO. 

BKL1*AKA 

BKS1*BKL1*CL1/CS1 
BKL2«BKLI*CL1/CL2 
C  IF  SIM«0.  THEN  A  CONPRESSIONAL  NAVE  IS  INCIOENT. 

SIN*|. 

CIM«CL1/<2.*PII 

CALL  COEF (BKS1.BKL1, 9KL2,SIN,CIN,RH0,GANAC) 

GAMAC*GAMAC*BKL1 

C  GAMAC  IS  THE  CONPRESSIONAL  SCATTERING  CROSS  SECTI0N*FRE1P*2. 
NRITE(6#55)  AKA, GAMAS , GAMAC 
55  FORMAT  f2XF7.2»2)5XE13 • A) ) 

GO  TO  1 
7  CONTINUE 
ENO 
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SUBROUTINE  C0EFIAKS1,  AKL1.  AKL2  ,S  I,  Cl, RHA, GAMMA* 

DIMENSION  *  J1  <400 )  f  AJ2(409I*  AJ3(60I>*  YH48»*t  Y2I400* 

THIS  SOLVES  THE  BOUNDARY  CONOITIONS  FOR  THE  AMPLITUDES  OF  THE 
SCATTERED  NAVES.  THE  AMPLITUDES  ARE  USED  TO  CALCULATE  THE 
SCATTERING  CROSS  SECTION. 

CALL  SPBES(AKS1«N1»AJ1) 

CALL  SPBESUKL1.N2,AJ2> 

CALL  SP8ES(AKL2»N3*AJ3) 

PI *3. 141592654 
IFIN1-N2)  10*11.11 

10  N*N2 
GO  TO  12 

11  N«N1 

12  IE<N-N3I  13*14*16 
13  NH*N3 

GO  TO  15 
16  NM*N 

C  NM  IS  THE  UPPER  LIMIT  OF  THE  SUM  IN  THE  SCATTERING  CROSS  SFCTION. 
15  CALL  SNNN(AKS1*NN,Y1) 

CALL  SNMN(AKL1«NM,Y2) 

GAMMAsO. 

00  100  M*l. NM 
F*FLOAT«HI-l. 

ZERO*C • 

011»E*AJ2IM)-ARLl*AJ2(M»l) 

012*ZERO-F*  (F»1.I*AJ1IM> 

RAT«F*U.-F»*AKSl**2/2. 

021»RAT*AJ2IM|-2.*AKL1*AJ2<M*1I 
122*F*IF*i.l*t  <F-l.»*AJltM»-AKSl*AJllH*l»* 
D3l«<F-l.»*AJ2(MI-AKLl*AJ2!H*ll 
ROT«l.-F**2*AKSt**2/2. 

032»R0T*AJHM*-AKS1*AJICN*II 

CIll«E*Y2IM)-A«Ll*r2IMU» 

CI12*F*«F*l.l*YltNI 
C13«F*AJ3IM)-AKL2*AJ3 IN*1> 

CI21*RAT*Y2M»-2.*AKLl*Y2<N»ll 

CI22*F*!F»l.)*ttF«l.»*UIH*-AKSl*Yl«N»ll» 

C23*RMA*AKSt**2*AJ3fNI/2. 

CI31«tF-t.»*V2M»-AKLl*Y2tM*lT 
CZ32*  R0T*YHM»-AKS1*YUN*1> 

Z1*(011*CI»012*SI> 

Z2«<021*CI*022*SI* 

Z3*(031*CI*032*SI ) 

A«N*022*C13*Z3*C23*032*ZI-C23*Z3*012-032*Z?‘C13 
AIM»C13*CI22*Z3*C23*CI32*Z1*C25*Z3*CI12-CI32*Z2*C13 
0RN«O31*Z2*C13»C23*Z3*Olt-C13*Z3*O21-C23*Zl*OSt 
9IN«C1S*Z2*CI31*C23*Z3*CI11-C13*Z3*CI21-C23*71*CI31 
ORH«C23* (012*031*CI12*CI31-0ll*032*CIll*CI32l  ♦CIS*  CD21*032*,CI?1* 
1CI32*C22*031*CI22*CI31 * 

0IM«C23*  f Cl 31*012*C  1 12*0  31 -C II 1*03 2-CI 32*0111  ♦CIS* I032*CI21*021* 
ICI32-CI22*031-022*CI31I 
P*ARN/0RM 
0«BRM/0RM 
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R*0IH/0RH 
S=AIH/ORH 
T*BIH/ORH 
0EN*1 . *R**2 

C  AH  AND  BH  ARE  THE  SQUARES  OF  THE  HAGNITUOE  OF  THE  SCATTEREO  HAVES. 
AH=IP**2»S**2I/0£N 
8H*IQ**2*T**21 /OEM 

GP*(2.  *F»1. 1  *4.*PI*  (A  *!/AKLl*F*  (Pf  i.»*BH/AKSi» 

GAHNA*GAMMA*GP 
190  CONTINUE 

GAHNA*8.6S6*GAHHA 

RETURN 

ENO 


> 

SUBROUTINE  SNHNI ZT,NN , Yl 
01 HENS  I ON  Y (A00) 

THIS  CALCULATES  THE  SPHERICAL  NEUHANN  FUNCTIONS  1 1 1 1  OF  ARG  ZT  AND 
UP  TO  OROER  NN.  THE  VALUES  ARE  CALCULATED  BY  UPHARO  RECURSION 
RELATIONS. 

ZERO*0. 

00  1  1*1, *00 

1  Y  C I) *0  . 

Y I II* ZERO-COS (ZTI /ZT 
Y(2l*ZERC-COS( ZTI /ZT*  *2-SINf ZT l/ZT 
NT*NN*2 
00  2  1*3, NT 
AJ*FL0AT <11-2. 

YUI*II2.*AJ*1,)/ZT1*Y  <1-1 1 -YU-21 
2  CONTINUE 
RETURN 
ENO 


SUBROUTINE  SP9ES<ZT,NT,AJI 
0IHENSI0N  RUIII,  RJIAOII,  AJ(AOO) 

THIS  CALCULATES  THE  SPHERICAL  BESSELS  FUNCTIONSfA Jf II I  OF  ARG 
ZT  UP  TO  OROER  NT. 

THIS  ROUTINE  IS  VALIO  FOR  ARGUICNTS  AS  LON  AS  O.BB. 

BT«ZT 

NDIH*A|| 

NZT*IFIX (ZTI 
00  lit  I *1 , NO  I H 
100  AU(II*0. 

THIS  SECTION  SETS  THE  UPPER  LI  HIT  OF  NT  FOR  THE  GIVEN  ARG  SUCH 
THAT  ALL  AJ#S  OF4  HIGHER  OROER  CAN  BE  SET  TO  ZERO. 

N*NZT+10 

00  1  I*N, 1850 , 5 
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J*I 

FI*FLOATtI> 

SECA*ZT/IFI».5) 

TACA=S0RTI1.-SECA**21 

C0SA=1./SECA 

SINA*SQRT(1.»C0SA**2) 

ALP*Al0GlCOSAfSQRTCC0SA**2-l.l » 

0ELT*EXP<  tFI*.5»*lTACA-ALPi  1/12.  •  «FI*.5)*SQRT!SINA1  I 
IF  C0ELT-1.E-27I  5,5,1 

1  CONTINUE 
GO  TO  12 

5  N* J 

FN*FLOAT«NI 

NP*N*1 

NT*N*2 

NU*N+1I 

EXPR*EXPtl.l 

00  r  I*NU, 1891,5 

K*I 

FI*FL0AT (II 
0*0.939299 

A1*2.*IFI-FN»1.I*0*AL0GCZT1 

A2*(2.*FI-2.*FN*1.»*0*ALOGIEXPR> 

A3*IFN*2. 51*0* ALOGIFN *2.1  ♦  <FN-1. 51 *0*ALOG!FN-2. 1 
A9*(FI*3.5»*O*AL0GIFI»3.) 

A5*IFI-.51*0*AL0GfFI-l.» 

E>m»«2M3-M*» 

IF(E*1I.)  9,5,7 
7  CONTINUE 

C  THIS  SECTION  CALCULATES  THE  BESSEL  FUNCTIONS. 

GO  TO  12 
9  NU*K 
J*NU 

R(NU*1I*0. 

00  2  1*1, NU 
K*I 

SJ«FLOAT(JI 

RUI«ZT/C1.*2.*SJ-ZT*R<J*11I 
IF(RCJ)-1.)  2,2,8 

2  J*J-1 

GO  TO  218 

6  IFU-2)  200,2)8,1% 

200  RJ (NU+1) *R(NU) 

RJtNU)*l. 

J*NU 

NUP*NU-1 
00  210  1*1, NUP 

J.J-I 

SJ*FLOAT(JI 

R  J(  J) *11 . *2.*SJI  /ZT*R J  <J*1I«RJ (J  +  2) 

210  CONTINUE 

ALPH«<RJ<l»-ZT*RJt2)l*COSI9mZT*SINI0T>*RJIll 
DO  220  1*1, NT 
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220  AJ<II*RJ<I»/ALPH 
333  CONTINUE 
RETURN 

14  RJ(JU>*RU> 

LAM2*«J*2 

IF(LAH2-NPI  15,44,44 

15  RJCJI *!• 

JsJ-1 

IC*K*1 

00  3  I*K,NU 
SJ*  FLOATUI 

RJ«J)*(l.*2.*SJ)/ZT*RJIJ*ll-RJtJ*2) 

3  J«J-1 

00  %  I*L AM2, NT 

4  Rja»*RJII-l»*RII-l) 

ALPH*(RJ«ll-ZT*RJI2»)*C0SnT|  ♦ZT*SINIBT»  •RJIll 
N0N*N0IH-1 
IF  CNT-NON)  16,16,17 
17  NT *N0M 

16  00  6  1*1, NT 

6  AJ(II*RJ(I)/ALPH 
105  CONTINUE 
RETURN 

44  WRITE (6,45) 

45  FORMAT I25H  LANBOA *2.GE.N»1  ^NO  XEQ) 

GO  TO  115 

12  MRITE (6, 131  ZT 

13  FORMAT <26H  ARG  9ES  FTN  TOO  LARGE  ZT*,E2I.8,5X14H  ENO  EXECUTION) 
STOP 

ENO 
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